Spin-polarized current oscillations in diluted magnetic semiconductor multiple 

quantum wells 
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We study the spin and charge dynamics of electrons in n-doped II-VI semiconductor multi- 
ple quantum wells when one or more quantum wells are doped with Mn. The interplay between 
strongly nonlinear inter-well charge transport and the large tunable spin-splitting induced by ex- 
change interactions with spin-polarized Mn ions produces interesting new spin dependent features. 
The tunneling current between quantum wells can be strongly spin polarized and, under certain 
conditions, can develop self sustained oscillations under a finite dc voltage. The spin polarization 
oscillates in both magnetic and nonmagnetic quantum wells and the time average in magnetic wells 
can differ from its zero-voltage value. Our numerical simulations demonstrate that the amplitude of 
the spin polarization oscillations depends on the distribution of magnetic wells within the sample. 
We discuss how the spin polarized current and the spin polarization of the quantum wells can be 
tailored experimentally. 

PACS numbers: 72.25.Dc, 72.25.Mk, 73.21.Cd, 75.50.Pp 
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I. INTRODUCTION 

Current information technology devices are based on 
the manipulation of electrical charge flow .-using electric 
fields. Recent technological developments^ that exploit 
magnetotransport effects in ferromagnetic metals have 
increased interest in exploring other ideas. Effects based 
on the electron spin degree of freedom can be used either 
to improve device functionality, or to create radically new 
devices that either implementnew processing algorithms, 
as in quantum computing, or thai-are based on new phys- 
ical principles, as in spin-FETSB The extreme sensitiv- 
ity to external magnetic fieldsu that has been exploited 
in ferromagnetic metals is ultimately a consequence of 
the collective behavior of many electronic spins that fol- 
lows from long-range ferromagnetic order. Ferromag- 
netic semiconductors are important because their mag- 
netic properties can in principle be engineered by doping 
or by adjusting gate voltages. Considerable progress has 
been recently made in manipulating the ferromagnetism 
that occurs in a number of common III-V compound 
semiconductors .when they are doped with the magnetic 
element MnffiE Although they are not normally fer- 
romagnetic, II-VI semiconductors doped with Mn, also 
respond stcongly to external magnetic fields and have 
propertiesBcm that are similar in many respects. The 
transport pr.aperties of Mn-based heterostructures have 
been studiedlj including miniband transport in strongly- 
coupled superlattices_doped with MnJiil 

In a recent paperE3 three of us analyzed the nonlin- 
ear transport properties of II-VI semiconductors with 
weakly-coupled highly n-doped quantum-wells in which 
one or more of the quantum wells is doped with Mn. The 
interplay between electronic spin, charge accumulation, 
and resonant interwell tunneling effects that takes place 
in these systems has been shown to alter the formation 



of stable electric field domains and to produce hysteretic 
behavior in both spin and charge distributions. Spin po- 
larization is induced in both magnetic and nonmagnetic 
quantum wells when driven by strong dc electric fields. 

In nonmagnetic highly n-doped semiconductor weakly 
coupled multiple-quantum- wells, it is well known that dc 
electrical transport is dominated by the formation of elec- 
tric field domains. This effect is reflected in the non- 
linear current-voltage characteristicsEjiiJ which present 
a sawtooth structure that arises from the interplay be- 
tween electron-electron interactions and resonant inter- 
subband tunneling between neighboring quantum wells. 
When the carrier density is below a critical value, how- 
ever, the formation of electric field domains is not stable 
and the non linear transport properties change drasti- 
cally in comparision with the highly n-doped systems. 
Intermediate n-doped systems have very rich behavior. 
The time-averaged current-voltage J—V characteristics 
presents flat plateaus. Real-time measurements show 
that within a plateau the cupent has an undamped oscil- 
latory dependence on time, il These self-sustained elec- 
tronic current oscillationsE3E3 come from the dynamics 
of the domain wall separating electric field domains and 
persist even at room temperature, making thesepdevices 
promising candidates for microwave generationllj with 
frequencies that extend from kHz to GHz. 

In this paper we explore new features that occur in 
weakly coupled quantum wells that support self- sustained 
oscillations when they are doped with magnetic impu- 
rities. In Section II we briefly review the model that 
we use for incoherent transport between quantum wells 
which can contain magnetic impurities. In Section III 
we summarize and discuss the results we have obtained 
by solving this model numerically for some representa- 
tive circumstances. Finally, in Section IV we present our 
conclusions. 
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II. MODEL 



Our theory of transport through diluted magnetic 
semiconductor multiple-quantum-wells is built on (i) 
a theory for the tunneling current between two spin- 
polarized two-dimensional electron gases (2DEGs); (ii) 
a continuity equation that accounts for relaxation of 
nonequilibrium spin populations; (iii) a relationship be- 
tween the up and down chemical potentials and their 
densities; (iv) a mean-field theory for the interaction be- 
tween 2DEG electrons and Mn spins; and (v) an expres- 
sion to take into account the Coulomb interaction be- 
tween charge accumulations in the quantum wells. 

In weakly coupled multiple-quantum- wells, it is a good 
approximation to treat tunneling of quasiparticles be- 
tween neighboring quantum wells by leading order per- 
turbation theory. We ignore interwell spin-flip processes, 
so that currents are carried between wells by the two spin 
subsystems in parallel. Accordingly, the current per spin 
a from the ith well to the (i + l)st well, Jfj + i, can be 
properly described by a Transfer Hamiltonian model00 
with a Lorentzian lifetime broadening of each quasi- 
particle's spectral density. Our Lorentzian halfwidth, 
7 = ?i/T sca tt, is energy independent, purely phenomeno- 
logical, and represents the combined effects of interface 
roughness, phonons, and impurity effects. Spin relax- 
ation procceses are not included in 7. The scattering 
times in quantum wells are shorter than any other time 
scale of the problem (r scatt — 0.1 — 1 ps); we therefore 
assume that the electrons in each well reestablish local 
equilibrium between succesive tunneling events and that 
their temperature is that of the lattice. The exchange in- 
teraction that couples s conduction band electron spins 
and d Mn local moments is ferromagneticQ and favors 
parallel alignment of the local moment S and band elec- 
tron spins s: H-f^t ~ Jsd Yli &i ' <%Ti) 1 where the sum 
is extended over the positions / of the magnetic impuri- 
ties and J s d is the exchange integral (which we take as a 
constant). 

When the mean-field and virtual crystal approxima- 
tions arc employed |l2 the effect of the exchange inter- 
action is to make the subband energies spin-dependent 
in those quantum wells that contain Mn ions: Ef = 



Ej - sA, where A = 2J sd N M nSB s (g^ B BS/k B T eS ) 
and s = +l/2(— 1/2) for a =TCl)- We must also take 
into account the microscopic processes that permit the 
quasiparticle system to bring its spin-subsystems into 
equilibrium within each .quantum well. Slow conduc- 
tion band spin relaxationt£l makes the effects we discuss 
stronger and is an important motivation for this study. 
Relaxation times in excess of 1 ns have been established 
experimentallycO in II- VI semiconductor quantum wells 
without Mn. In II- VI diluted magnetic semiconductor 
quantum wells these times are reducecL-to tens of picosec- 
onds (but are still larger than r sca tt).cll Fermi's Golden 
rule leads to the following equations for the spin relax- 
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FIG. 1: J—V curves for a 50-well n-doped 5-nm ZnSe/10- 
nm Zno.8Cdo.2Se system with fractional MnSe monolayers at 
the 1st and 50th quantum wells. Well doping is (a) N w = 



1 x 10 1 



and (b) N w = 9 x 10 1 



Contact doping 



is (a) iV c = 1.1 x 10 12 cm" 2 and (b) N c = 9.9 x 10 10 cm" 2 . 
Case (a) corresponds to a highly n-doped sample. In that 
case, the current is stationary and electric field domain forma- 
tion takes place. It is reflected in the sawtooth-like structure 
that presents J versus V. Case (b) corresponds to an interme- 
diate n-doped sample. In that case, the current presents self- 
sustained oscillations and the time-averaged current shown in 
the figure above shows a flat plateau behaviour. 



ation rate equation within each well: 



dn\ 
~~dt 



T s f 
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(1) 



where W is the spin opposite to er, r s f is the spin-scattering 
time, vq is the two-dimensional density-of-states per spin, 
and nf is the spin a charge density in the ith-well. For 
A > /xj - E\ x , Eq. (0) must be modified!! 
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For large enough A, Eq. (g) leads to an equilibrium state 
with full spin polarization. The continuity equations that 
describe the time evolution of the partial densities in each 
quantum well include both spin-relaxation and transport 
currents: 



dt 



U i—X,i d i,i+l 



1/0 i=l,...,N (4) 



for the case fij — Ej 1 > A (N is the number of quantum 
wells). Otherwise, Eqs. (||||) must replace the second 
term on the right-hand side of Eq. (0) . 
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III. RESULTS 



In order to obtain electric field domain physics it 
is necessary to account for electron-electron interac- 
tions among conduction band electrons at— the mean- 
field level by solving the Poisson equation.^ Including 
a displacement current contribution, the total current 
density J(t) traversing the sample at time t is J(t) = 
2^W ^M+iW which is independent of i and where 



Ji 



4 l+ l (*)" 



Ji i+ i(t), e is the sample permittiv- 



ity, d is the multiple-quantum well period and Vi the volt- 
age drop between wells i and i + We model a N — 50 
n-doped ZnSe/Zno.gCdo^Se multiple-quantum-well sys- 
temEl with well and barrier widths 10 nm and 5 nm, re- 
spectively. Mn impurities are placed in the 1st and 50th 
quantum wells. We have chosen the spin relaxation time 
within the magnetic and nonmagnetic quantum wells to 
be T s f — 10 -11 s and r B f = 10~ 9 s, respectively. The 
quantum wells and the contacts are n-doped with N w 
and N c = kN w respectively. We take re = 1.1. Figure fil 
shows the J-V characteristics for (a) N w = 1 x 10 12 cm 
and (b) N w = 9 x 10 10 cm" 2 . The spin splitting is A = 2 
meV in both cases. With these parameters, the magnetic 
quantum wells (if taken as isolated) are partially polar- 
ized: P = 7% [case (a)] and P = 75% [case (b)], where 
Pi = (nj —n\)/rii is the spin polarization of the zth-wcll. 
As expected, Fig. 0(a) exhibits a series of sharp discon- 
tinuities in the negative differential conductance region 
which are linked to the formation of static electric field 
domains. Along each branch, charge accumulates at the 
domain wall, forming a monopole that jumps discretely 
toward the emitter as V increases. When N w is lowered 
[see Fig. 0(b)], the branches are substituted by a plateau 
on which current oscillates periodically with time. We 
observe undamped self-oscillations of the current in the 
range V = 0.2-1.5 V. 

In Fig. |we plot the total current (a), and the spin-up 
(b) and spin-down (c) currents toward the collector as a 
function of time for V = 0.5 V for the sample described 
in Fig. |l](b), i.e., for an intermediated n-doped sample: 
N w = 9 x 10 10 and N c = 9.9 x 10 10 which presents self 
sustained oscillations in the absence of magnetic quan- 
tum wells. The current oscillations have a period of the 
order of 5 /is, i.e., with a frequency in the kHz range. We 
observe an irregular shape for the current amplitude in all 
the three cases. The current spikes that are superimposed 
on the main periodic structure reflect discrete jumps of a 
domain wall separating electric field domains from well to 
well.t£l The current oscillations emerge from the dynam- 
ics of the domain wall which usually consists of a charge 
accumulation layer (monopole) spread over one or two 
quantum wells. In nonmagnetic multiple-quantum-wells 
the dynamics of the monopole has been theoretically de- 
scribed and experimentally observed.Lj Lowering the con- 
tact charge doping, the domain walls lead to traveling 
dipolesE3 that consist of one accumulation and one deple- 
tion charge layer located in separated quantum wells and 
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FIG. 2: (a) Total time-dependent current (tunnel- 
ing+displacement), (b) spin-up, and (c) spin-down time- 
dependent current at a fixed dc bias voltage Vd c = 0.5 V 
for the superlattice described in Fig. [j](b). The current oscil- 
lations present a flat region and overimposed spikes (see the 
text below). Comparison of (b) and (c) indicates that the 
current towards the collector is partially spin-up polarized. 
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FIG. 3: Local distribution of the electron density through 
the sample at fixed dc bias voltage: Vdc = 0.5 V for the 
different times marked in Fig. 0(a). At the flat region of the 
oscillation in Fig. 0(a) the domain wall consists of one charge 
accumulation layer (monopole) which propagates just within 
a small part of the superlattice. With increasing time (ts — 
tio) the domain wall configuration changes and it becomes 
a dipole (one charge accumulation and one charge depletion 
layer) which runs through half of the superlattice. 



have distinct dynamics. (The dipole charge front wave is 
similar to the one responsible of the Gunn effect in bulk 
semiconductors). In semiconductor multiple-quantum- 
wells, monopole. and dipole domain walls can coexist at 
a fixed dc biasE^I under certain conditions. 

In Fig. 0(a) we plot the charge density quantum well 
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FIG. 4: Spin polarization, P, in the quantum well closest to 
the collector as a function of t for Vdc = 0.5 V. The frac- 
tional polarization of the isolated quantum well is P = 75%. 
Within the superlattice, in the strongly nonlinear regime, the 
polarization oscillates and reaches, for a small time window 
of the period, full spin-up polarization. 
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FIG. 5: Spin polarization of the magnetic quantum well, lo- 
cated in the 10th (a), 20th (b), 30th (c), and 40th (d) posi- 
tion at a fixed dc bias, Vdc = 0.5 V. In all four cases a second 
magnetic quantum well is placed at the location closest to 
the emitter. The polarization corresponding to the isolated 
magnetic well is P — 75%. The depletion charge of the dipole 
domain wall is centered around the 40th well. This depletion 
is responsible for the enlargement of the spin polarization and 
can lead to full polarization [see (d)]. 



distribution at the times ti-t$ marked in Fig g. We ob- 
serve that for the case illustrated the domain wall is a 
monopole which travels through only a part of the struc- 
ture. In Fig. ||(b) the charge density at times ts-tio is 
presented. These plots illustrate an interesting feature of 
our results: the domain wall is now a dipole with both 
charge accumulation and depletion. During each oscilla- 
tion period the domain wall undergoes a transition from 
a monopole to a dipole one and back. An exciting con- 
sequence of this behavior is shown if Fig. [| where the 
polarization in the magnetic quantum well is plotted as 
a function of t. We observe that the polarization reaches 
three different values during one oscillation: during the 
intervals of dipole propagation the polarization abruptly 
drops and increases up to a practically fully polariza- 
tion configuration. When the domain wall becomes a 
monopole the polarization remains practically constant 
up to the time where the dipole is formed. The con- 
stant value of the polarization is close to the value of the 
polarization of the isolated quantum well. 

The previous observation can be explained by looking 
at how the spin polarization in the Mn-doped quantum 
well influences the tunneling probability to neighboring 
wells and by the strongly nonlinear transport in the neg- 
ative differential conductance region. The chemical po- 
tentials for spin up and down electrons within a magnetic 
quantum well do not take on the same value in the sta- 
tionary limit as they would do in an isolated quantum 
well because of the finite spin-flip and tunneling rates 
and because self-consistency changes the charge density 
in the quantum wells. 

The magnetic quantum well spin-polarization in the 



nonequilibrium case may be related to the instantaneous 
local chemical potentials for majority and minority spins: 

P = 1 , (5) 

where S^, m — fi^ — This expression is obtained as- 
suming the linear non-interacting 2DEG relationship be- 
tween density and chemical potential, so that correlation 
effects are not included:E 3 nf = - Ef x ). The above 

expression for the polarization P is also valid for non 
magnetic quantum wells. In this case, the spin splitting 
is very small and it is basically the second term of Eq. (H) 
that is responsible for the induced finite polarization. 

As we can see from Eq. (|j), increasing the den- 
sity reduces the relative polarization. It follows that a 
monopole domain wall accumulation layer in the mag- 
netic well will reduce P compared to the equilibrium case. 
This reduction was observed for multiple-quantuipiwells 
with stable electric field domains and stationaryt3 cur- 
rents. In contrast, in the self-sustained current oscilla- 
tions the magnetic well polarization can increase with 
respect to the isolated case, even reaching full spin- 
polarization. Comparing Figs. || and ^, we note that the 
high polarization configuration state occurs at t = tg, 
when the front wave is a dipole with one depletion and 
one accumulation layer (domain walls) and the depletion 
region is located at the magnetically doped quantum well. 
The large polarization occurs in part because of the local 
instantaneous decrease in the total density of electrons in 
the magnetic well. The frequency, amplitude and shape 
of the current oscillations depend on the dynamics of the 



5 



polarized charge, which itself depends on the distribu- 
tion of magnetic quantum wells in the sample. Then 
tailoring the sample configuration should enable control 
of not only the frequency and amplitude of the current 
oscillations, but also of the spin polarization and its time 
dependence. To illustrate this we present in Fig. || re- 
sults for the case in which two magnetic quantum wells 
are placed in the system, one adjacent to the emitter and 
another at a variable position. We observe that when the 
magnetic quantum well is close to the position of the do- 
main wall the amplitude of the polarization oscillations 
increases [Fig. ||(d)]. Because the polarization depends 
on the inverse of the charge density, whenever the do- 
main wall enters and leaves the magnetic quantum well, 
the amplitude of the magnetization oscillations increases 
and P oscillates between 20% and 100%. As the magnetic 
quantum well comes closer to the collector, the polariza- 
tion exhibits two flat regions within an oscillation. Since 
the period of these electric field domain oscillations is of 
the order of 5 /xs, time resolved measurements of PL in 
illuminated multiple quantum well samples should allow 
this polarization change to be observable. 

IV. CONCLUSION 

In closing, we have studied the spin dynamics in 
multiplc-quantum-well systems doped with magnetic im- 



purities. We predict time-dependent periodic oscillations 
of the spin polarized current and of the spin polariza- 
tion in both magnetic and nonmagnetic quantum wells. 
The interplay between nonlinearity of the current-voltage 
relationship and the exchange interaction produces new 
spin dependent features which are sensitive to spin re- 
laxation times and to the equilibrium exchange fields in 
the quantum wells which can be tailored by adjusting ex- 
ternal field and temperature. These new spin dependent 
features can potentially be exploited for device applica- 
tions. In particular, our results for the oscillating spin 
polarized currents at the collector suggest that these sys- 
tems could be designed as spin-polarized current injection 
oscillators. 
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